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INTRODUCTION 


Hybopsis  harperi  Is  a  rheotropic  cyprlnld  minnow  found  in 
Florida  in  two  extremely  different  situations  and  in  several  variants  of 
each.   Although  spring-fed  runs  and  rivers  constitute  its  main  habitat, 
H.  harperi  is  also  found  in  small,  steep- sided  sink  holes  in  the  Ocala 
limestone  of  North  Central  Florida.    Hellier  (1957),  in  his  study  of  the 
Santa  Fe  River,  pointed  out  that  It  Is  the  most  abxuidant  fish  In  the 
spring-fed  stretches  of  the' river.   Marshall  (1941,  1947)  discussed  the 
seemingly  strong  ecologlc  dlvergeac«  of  the  two  habitats  occupied. 

Marshall  considered  the  absence  of  Gambusla  from  most  sink 
holes  as  evidence  for  dispersal  of  Hybopsis  through  a  subterranean 
water  system.    The  presence  of  a  large  population  of       harperi  in  a 
pond  In  a  lime  rock  mine  near  Halle,  Florida,  further  suggests  this 
means  of  dispersal.    Within  a  year  of  Us  filling  with  water,  this  pond 
contained  Hybopsis  but  no  Gambusla,  even  though  it  was  undoubtedly 
suited  for  occupation  by  Gambusla.    The  ecological  suitability  of  the 
habitat  is  Implied  by  the  fact  that  a  local  bass  fisherman  had  succeeded 
in  Introducing  Xlphophorus  variatus  into  the  pond  In  order  to  have  a 
ready  suppJ"  of  bait.    Xlphophorus,  a  poeciliid  not  native  to  Florida, 
appeared  to  be  occupying  the  same  niche  that  Gambusla  does  In  sink 


holes  Into  ^Ich  It  has  been  Introduced. 

If  dispersal  of  H.  harperi  occurs  through  the  underground  system, 
there  is  certainly  some  degree  of  isolation  among  the  subpopulations  in 
the  sink  holes,  and  also  between  these  and  the  ones  in  streams.  In 
fact,  it  is  conceivable  that  in  some  cases  there  may  be  a  greater  degree 
of  Isolation  between  neighboring  sink  holes  than  between  these  and  some 
distant  stream* 

Therefore,  unless  we  are  to  postulate  some  directed  migration 
among  the  various  populations ,  H.  harperi  in  northern  Florida  is  evi« 
dently  a  large  popxilation  composed  of  semi>lsolated  subpopulations. 
With  this  population  strixcture  we  woidd  expect,  in  the  presence  of 
selection  pressure,  to  find  adaptive  changes  in  the  cave  and  sink  hole 
populations.   Moreover,  under  Wright's  (1931)  hypothesis,  these  changes 
should  occur  rather  rapidly.   In  the  absence  of  selection  pressure,  ran- 
dom changes  In  gene  frequencies  would  be  escpected  to  manifest  th em- 
salves  la  altaxations  of  morphometrlc  or  merlstic  characters.  Also, 
with  the  degree  of  homozygosity  which  results  from  Inbreeding  In 
isolated  subpopulations,  there  might  be  an  increase  in  nongenetic  vari- 
ability (Lerner.  1954). 

Thus,  there  are  two  factors  which  could  be  expected,  a  priori,  to 
contribute  to  the  formation  of  subspecifically  distinct  populations. 
First,  one  can  offer  H.  harperi  as  a  model  for  the  situation  in  which 


Wright  (1931)  suggested  that  evolution  would  proceed  most  rapidly; 
i. «.  *  a  large  population  broken  up  into  small,  almost  (but  not  com- 
pletely) isolated  units.    Secondly,  if  there  is  no  migration  at  all  be- 
tween the  subpopulations ,  one  would  expect  genetic  drift  to  have 
allowed  some  nonadaptive  differences  to  develop* 

Under  the  first  of  these  hypotheses,  one  can  even  postulate  a 
direction  for  this  presumptively  rapid  evolution.    Sink  holes  such  as 
those  included  in  this  study  form  a  graded  series  from  fairly  well- 
lighted  ones  to  completely  dark  caves.    Thus  the  population  structure 
ef  Hybopsis  could  assist  in  the  bridging  of  the  gap  between  the 
"adaptive  peak"  corresponding  to  the  stream  habitat  and  that  corre- 
sponding to  the  cave  or  sink  hole  habitat.    No  blind  or  unpigmented 
Hybopsis,  which  might  be  expected  if  evolution  had  proceeded  far 
enough  in  the  direction  suggested,  have  been  found. 

In  addition  to  the  nonadaptive  divergence  of  populations  caused  by 
random  genetic  drift,  "accidents  of  sampling"  might  also  be  expected 
in  the  colonisation  of  sink  holes,  or  in  the  re-establishing  of  popula- 
tions from  survivors  of  periodic  reductions  during  times  of  drought. 

Regardless  of  the  causes  of  variation  in  Hybopsis  harperi,  there 
would  seem  to  be  grounds  to  expect  that  a  study  of  this  variation 
would  reveal  differences  sufficient  to  justify  separation  of  at  least 
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two  subspecies.   Hubbs  and  Crowe  (1956)  did  this,  describing  the  sink 
hole  form  as  a  distinct  subspecies,  Hybopsis  harperi  subterranea. 
This  fact  and  the  pattern  of  local  distribution  of  Hybopsis  in  northern 
Florida  suggest  that  a  further  study  of  variation  in  this  species  might 
shed  light  on  the  effect  of  poptdation  structure  on  subspeciation  and  on 
the  course  of  microevolution  generally. 

Conventional  practices  with  respect  to  the  naming  of  subspecies 
have  been  subjected  to  intensive  examination  in  recent  years.  Be- 
ginning with  a  paper  by  Wilson  and  Brown  (1953),  a  series  of  some- 
what polemic  essays  has  appeared  in  Systematic  Zoology,  with  such 
proposals  as  that  trinomens  be  stripped  of  both  priority  and  italics 
(Gosline,  1954)  or  that  all  subspecies  be  completely  abolished  (Brown 
and  Wilson,  1955).    There  has  even  been  a  recommendation  that  an 
"objective"  subspecies  definition,  based  solely  on  geographic  con- 
siderations, be  adopted  (Edwards,  1956),  but  Borgmeier  (1957) 
favored  a  morphologically  conceived  subspecies.    Hubbell  (1954), 
summarizing  a  symposium  on  the  subject,  concluded  that  in  general 
Infraspecific  variation  is  best  treated  by  description,  graphic  repre- 
sentation, and  nontechnical  names,  and  not  by  formal  names  subject  to 
the  rules  of  nomenclature. 

Except  for  some  purely  legalistic  questions  of  nomenclature,  cri- 
ticism of  the  current  subspecies  concept  seems  to  center  arovmd  three 


main  Issues,  as  follows:  (1)  the  criteria  for  distinguishing  subspecies 
are  usually  quantitative  characters  v^ich  do  not  exhibit  complete  dis- 
continuity; (2)  analysts  of  several  characters  fretjuently  shows  that 
there  is  nonco&cordance  in  the  variation  of  these  characters;  and 
(3)  the  adoption  of  a  "percent  rule,"  ti^ich  is  made  necessary  by  the 
lack  of  discontinuity,  is  altogether  arbitrary.  ' 

Two  recent  publications,  neither  of  them  in  the  mainstream  of  the 
subspecies  controversy,  have  dealt  with  the  statistical  aspects  of  the 
problem.   Plmentel  (19S8)  suggested  the  study  of  many  characters  and 
material  from  many  tocallttes  before  reaching  a  conetusloa  concerning 
the  probable  existence  of  subspecies.  He  also  outlined  a  sound  statis- 
tical procedure  for  such  a  study.   Marr  (1957)  defined  several  levels 
of  population  units  In  fishes  and  proposed  methods  for  studying  the 
subpopulatlon  s . 

The  present  stvuly  stands  somewhere  between  these  two  papers  In 
its  approach  to  the  problem  of  subspeclatlon  In  Hybopsis.  Plmentel, 
who  was  concerned  primarily  with  morphological  data,  was  especially 
interested  in  the  assured  randomness  of  samples.  Although  he  en- 
couraged Inclusion  of  many  characters,  his  methods  were  restricted 
to  analysis  of  each  character  separately.   Marr,  on  the  other  hand, 
extended  his  techniques  of  analysis  to  include  biochemical  and  tagging 
data  and  genetic  studies  wherever  they  are  possible.   He  also 


suggested  the  use  of  multivariate  analysis  la  the  morphometrlc  phase 
of  an  Investigation,  and  Royce  (1957) »  In  the  same  symposium,  de- 
veloped this  technique.  , ;  , 

My  study  has  utilized  some  specimens  from  museum  collections 
and  others  collected  In  much  the  same  way  as  those  that  furnish  the 
material  of  the  usual  Ichthyologlcal  Investigation.    Ilierefore^  In 
contrast  to  Plmentel's  suggestion*  randomness  of  samples  Is  not 
assured.   Even  though  the  populations  of  variables  are  thus  non- 
randomly  sampled,  there  Is  no  reason  to  believe  that  the  departvires 
from  randomness  are  the  same  for  all  variables;  so  the  synthesis  of  a 
composite  multivariate  statistic,  such  as  Mahalanobls'  D^,  should  be 
less  affected  by  the  nonrandomness  than  single  variable  statistics. 
My  method  of  evaluating  varlatlcm  In  subpopulatlons  of  Hybopsls 
harperl  has  been  to  use  tiie  analysis  of  variance  technique,  as 
suggested  by  Plmentel«  on  single  characters  and  then  to  extend  the 
analysis  to  multivariate  methods,  as  suggested  by  Marr  and  by  Royce. 


,  !    MATERIALS  AND  METHODS     -  -  • 
:.    Sample  a 

The  data  compiled  for  ^Is  sttuiy  conslat  of  counts  and  measurements 
from  a  total  of  635  specimens  of  Hybopsis  harperi ,  all  collected  Ln 
North  Central  Florida  (Figure  1).    Some  of  the  specimens  were  in  the 
University  of  Florida -Florida  State  Museum  collections;  the  rest  were 
collected  by  the  author.   For  the  purpose  of  analysis  it  will  be  con- 
venient to  recognise  three  collection -lots  in  the  mjiterial:  a  Special 
Series*  a  Santa  Fe  Series*  and  a  Main  Series.   In  the  following 
gazetteer  of  localities  for  these  collections,  the  abbreviation  by  which 
each  will  hereinafter  be  designatod  is  given  first.   The  first  two  letters 
of  the  abbreviation  indicate  the  location,  as  shown  in  Figure  1;  the 
third  letter  indicates  whether  the  habitat  is  standing  water  or  running 
water,  and  a  number  designates  one  of  several  samples  collected  at 
the  same  place. 

The  Special  Series  comprises  a  random  sample  of  100  specimens 
from  each  of  the  following  three  collections: 

SFK-I:    A  collection  (UF-8425)  made  by  Hellier  in  the 
Santa  Fe  River  near  Duncan's  Ijanding  on 
November  17,  1956.  ' 


HPS:      A  collection  made  from  a  lime  rock  mine  near 

Haile  on  January  29,  1959.    Thi«  pool  became  pop- 
\ilated  by  Hybop«i»  within  a  year  after  the  mining 
operation  reached  the  water  table,  and  the  environ- 
ment it  quite  unlike  the  stream  habitat  or  the  sink 
hol«  habitat. 

ZMS:     A  collection  from  Zamia  Sink,  a  sink  hole  in 
Alachua  County,  made  on  December  29,  1958. 
The  Santa  Fe  Series  is  made  up  of  five  more  samples  from  the 
same  section  of  the  Santa  Fe  River  as  SFR-1.    These  collections,  also 
made  by  Hellier  during  the  course  of  his  study  of  the  Santa  Fe,  are  as 
follows:  ^ 

SFR'Z:  A  collection  made  February  25,  1956,  by  means  of 
a  common  sense  seine  and  a  ten-foot  bag  seine. 

SFR-3:  A  collection  made  July  3,  1955,  with  a  comtmon 
sense  seine  and  a  ten-foot  bag  seine. 

SFR-4:  A  collection  made  July  7,  1956,  with  a  fifty-foot 
bag  seine  and  a  conomon  sense  seine. 

SFR-5:  A  collection  made  on  July  22,  1956,  with  a  fifty- 
foot  bag  seine. 

SFR-6:  A  collection  made  September  14,  1956,  using  a 
fifteen- foot  bag  seine. 
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Collection  SFR-3  contained  only  25  specimen*,  all  of  which  were 
included  in  the  sample .   Each  of  the  others  was  a  random,  sample  of 
30  specimens  from  one  of  the  above  collections. 

The  Main  Series  is  made  up  of  samples  from  a  variety  of  stream 
and  sink  hole  environments  in  North  Central  Florida.    Each  is  a  ran- 
dom sample  of  30  specimens  from  one  of  the  collections  indicated 
below:  j  ■  ■  ■  - .  ■ 

JRS-1:    A  collection  made  on  July  19,  1959,  from  Jerome 
Sink,  a  sink  hole  in  Alachtia  County.   This  is  the 
type  locality  of  Hybopsis  harperi  subterranea 
Hubbs  and  Crowe. 
JRS-2:    A  collection  from  Jerome  Sink  which  was  made 

May  30,  1955. 
CWS:      A  collection  from  Cow  Sink  which  was  made 

January  29,  1959.   Cow  Sink  is  located  about  200 
yards  west  of  Jerome  Sink. 
BLR:      A  collection  from  Blue  Springs  in  Gilchrist 
County  on  June  11,  1959.   Blue  Springs  Is  a 
tributary  of  the  Santa  Fe  River. 
SGR:       A  collection  made  from  Silver  Glen  Springs  on 
September  20,  1951  (UF-6307). 
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SLR:     A  collection  from  Salt  Springs  on  November  9,  1949 
(UF-6303). 

•  JUR:     A  collection  from  Juniper  Springs  made  by  Dr. 

JohnD.  Kilby  on  November  1,  1958. 
Counts  and  Measurements  ^ 
The  values  for  five  meristic  characters  were  determined  for  the 
specimens  in  the  Special  Series  and  in  aU  samples  of  the  Main  Series 
except  SGR  and  SLR.    These  characters  and  their  abbreviations  are  as 
follows:  dorsal  fin  rays  (D),  anal  fin  rays  (A),  scales  in  the  lateral 
line  (LL),  scales  above  the  lateral  line  (SA),  scales  below  the  lateral 
line  (SB).   An  analysis  of  variance  for  the  meristic  characters  is  pre- 
sented below.    Counts  were  discontinued  in  the  later  stages  of  the 
study  after  a  preliminary  analysis  of  the  Special  Series  showed  no 
reliable  differences  in  the  meristic  characters  in  these  three  samplss. 
Values  for  thirteen  morphometric  characters  were  obtained  for  each 
specimen.    The  meristic  counts  were  taken  according  to  Hubbs  and 
Lagler  (1941),  and  wherever  possible  the  measurements  were  taken  be- 
tween the  reference  points  suggested  by  these  authors. 
The  morphometric  characters  are  as  follows: 
AG-------Angle  of  gape 


L  Standard  Icogth 

DB  Depth  of  body 

HLi- ---- Head  length 

SL  Length  of  snout 

E  Width  of  eye 

LD  Length  to  dorsal  fin 

PD----- Depth  of  caudal  peduncle 
HW  Width  of  head 

HL/SL --Length  of  head/length  of  snout 

HL/E-— Length  of  head/width  of  eye 

HL/HW-  Head  length/head  width 

PD/DB- -Depth  of  peduncle/depth  of  body. 
All  the  linear  measurements  except  L,  DB,  and  LD  were  mad* 
with  the  aid  of  a  camera  luclda.   The  image  of  the  fish  was  projected 
onto  a  sheet  of  white  paper  on  which  index  lines  had  been  previously 
Inscribed,  and  the  reference  points  were  marked  on  the  index  lines. 
The  measurements  were  then  taken  from  the  sheet  with  dividers. 
The  magnification  of  the  projected  image  was  14.  2  diameters.  This 
method  was  not  applicable  for  L.  DB,  or  LD  because  the  reference 
points  for  these  measurements  could  not  be  included  in  one  field.  An 
opaque  projector,  magnifying  1.57  diameters,  was  used  to  project 
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these  reference  pointe  onto  another  index  line  on  the  same  sheet. 
These  points  were  also  marked  and  the  measurements  taken  with 
dividers  and  the  same  scale. 

For  the  purpose  of  this  statistical  study,  the  actual  values  of  the 
measurements  were  not  determined;  the  computations  were  per- 
formed on  measurements  in  fiftieths  of  an  Inch*  taken  directly  from 
the  data  sheets.   To  convert  these  to  acttial  metric  values  In 
millimeters,  one  may  multiply  JL,  DB,  and  JLD  by  0.3236,  or  the 
other  measurements  by  0.0358.    Obviously  the  statistical  results 
would  not  change  with  this  transformation.   In  order  to  facilitate 
comparisons,  Table  7  expresses  the  means  and  standard  deviations  of 
all  samples  In  millimeters. 

The  reference  points  for  the  angle  of  gape  (AG)  were  also  marked 
by  projection  through  the  camera  luclda.    The  angle  measured  was  that 
formed  by  the  "horlaontal"  index  line  and  a  line  drawn  between  the 
anteriormost  and  posterlormost  extent  of  the  gape.   The  difficulty  of 
making  a  true  horizontal  superlmposltlon  of  the  Index  line  produced 
rather  large  variations  in  this  measurement.     The  values  for  AG 
reported  here  probably  can  not  be  compared  with  those  reported  by 
Hubbs  and  Crowe  (1956)  because  of  differences  In  the  method  of 
measuring  the  angle. 


"  Computationa 
Computations  of  the  sums  of  sqiiares,  stuns  of  products,  means » 
and  standard  deviations  for  each  sample  were  performed  on  the  IBM 
650  computer  at  the  University  of  Florida  Statistical  LAboratory, 
These  were  summarised  by  use  of  the  407  accounting  machine  at  the 
IBM  Supplies  Oivisioa  Card  Plant  at  Sherman,  Texas.  Inversion  of 
the  varlance-covarlance  matrix  was  done  by  the  Statistical  Laboratory, 
as  were  the  matrix  multiplications  necessary  for  the  computation  of 


SINGLE  FACTOR  ANALYSIS 
Analysis  of  Variance  . 

Since  multivariate  statistical  analysis  re^viires  computations  that 
are  laborious  even  on  modem  electric  desk  calculators,  single  factor 
analyses  should  be  used  ii  they  will  allow  valid  conclusions  to  be  drawn. 
Therefore,  the  analysis  of  variance  technique  was  used  here  In  an 
attempt  to  identify  any  of  the  characters  which  might  allow  discrim- 
ination among  the  subpopulatlons  on  the  basis  of  single  factors.  In 
addition  to  the  ustial  analysis  of  variation  between  samples  generally, 
I  have  also  been  able  to  make  orthogonal  comparisons  in  a  search  for 
the  source  of  variations  between  samples.   The  one  odd-sized  sample 
in  the  Santa  Fe  Series  prevents  complete  orthogonality  of  any 
meaningful  comparisons.    Thus,  the  fourth  comparison  Is  merely  a 
measure  of  the  variation  remaining  after  the  other  comparisons  have 
been  made.    That  a  standard  sample  size  shoxild  be  adopted  In 
variational  Investigations  Is  Indicated  by  the  fact  that  much  of  the 
significant  variation  was  In  that  remainder. 

Although  merlstlc,  scalar,  linear,  and  angular  characters  may  be 
Included  in  the  same  multivariate  analysis  (Olsen  and  Miller,  1958), 
preliminary  analyses  of  tiie  Special  Series  and  the  differences  In  units 
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led  to  a  decision  to  restrict  the  analysis  o£  the  meristic  characters  to 
the  single  factor  analysis  of  variance,  including  orthogonal  sub- 
division of  the  between- sample  siuns  of  squares. 

One  of  the  assumptions  in  the  analysis  of  variance  Is  that  the 
variances  In  the  different  samples  are  homogeneous.  Departures 
from  this  assumption  will  result  in  a  loss  of  sensitivity  in  the 
analysis  (Cochran,  1947).   Even  though  this  technique  was  to  be  used 
primarily  for  Identifying  characters  showing  variation,  Bartlett's 
test  of  the  homogeneity  of  variances  (Ostle,  1954)  was  performed.  In 
all  three  series,  the  hypothesis  of  homogeneity  of  the  variances  was 
rejected  for  most  of  the  characters  (Table  8,  Appendix).   Li  (1957), 
Snedecor  (1956),  and  Cochran  and  Cox  (1957)  suggest  that  moderate 
departures  from  homogeneity  are  not  serious  in  analysis  of 
biological  material,  but  they  all  caution  the  worker  to  be  careful  In 
Interpreting  the  significance  level  of  the  differences.   In  this 
cozmectlon.  It  Is  Interesting  to  note  that,  of  161  F  values  which 
exceeded  the  5  percent  tabular  vadue  in  these  variance  analyses,  ISO 
also  exceeded  the  0.1  percent  level.    One  may  then  cautiously  s\iggest 
that  these  F  values  indicate  real  differences  in  tiie  samples  on  which 
they  are  based.    Throughout  this  report  the  values  which  exceed  the 
5  percent,  1  percent,  and  O.l  percent  values  for  T   will  be  Indicated 
respectively  by  single,  double,  and  triple  asterisks. 
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TABLE  1 


INTERPRETATION  OF  THE  ORTHOGONAL  COMPARISONS 


C<»npari«on  Series 


A  significant  value  indicates  a 
difference  between  the: 


C, 


\"  Santa  Fe 


C4 


Special  sample  from  running  water  and  the  two 

->y*if^        '..i:    samples  from  still  water. 

Special  two  samples  from  standing  water. 

samples  collected  with  a  fifty- foot  bag 
seine  and  those  taken  with  shorter 
..  .,  seines.   

Santa  Fe        two  samples  collected  in  the  same  area 
during  the  same  month.  . 

Santa  Fe        samples  collected  in  February  and  the  one 
taken  in  September. 

Santa  Fe        samples  in  this  series  which  is  not 

accounted  for  by  the  other  three  com- 
-  .  parisons. 

Main  samples  from  springs  and  those  from  sink 

holes.  :  1  , 


^10 


Main 


Main 


Main 


samples  from  two  closely  adjacent  sink 
holes,  Jerome  Sink  and  Cow  Sink. 

two  samples  taken  from  Jerome  Sink  at 
different  times. 

sample  from  Blue  Springs  and  the  three 
from  Marion  County  springs.    (For  the 
meristic  characters,  this  comparison  is 
only  between  Blue  Springs  and  Juniper 
Springs.) 
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TABLE  1- Continued 


Comparison  Seriea 


A  aignificant  F  value  indicates  a 
difference  between  the: 


Main 


sample  from  Juniper  Springs  and  the  other 
two  from  Marion  County. 


Main 


sample  from  Salt  Springs  and  that  from 
Silver  Glen  Springs.  


In  the  Special  Series,  two  single-degree-of-freedom  orthogonal 


comparisons  were  possible.   Sums  of  squares  for  these  and  all  other 
comparisons  were  computed  using  the  method  of  coefficients  (Ostle, 
1954;  p.  267).    The  first  Special  Series  comparison  (indicated  by  C^) 
tests  the  differences  in  the  characters  which  are  related  to  the 
differences  in  standing-  or  running-water  habitats.    The  second  com- 
parison (C2)  compares  the  two  standing -water  samples,  from  Zamia 
Sink  and  the  lime  rock  pit  at  Haile.   A  similar,  but  more  extensive, 
analysis  could  be  made  for  the  Main  Series.    In  this  series,  Cy  ' 
measures  differences  between  springs  and  sink  holes;  Cg  compares 
two  sink  holes;  C9  compares  samples  from  the  same  sink  hole  in 
different  years;  and  Cio  compares  two  springs.    The  notation  for  all 
orthogonal  comparisons  is  consistent  throughout  this  report,  and  the 
coefficients  for  these  comparisons  are  presented  in  the  Appendix, 
Table  9.    Table  1  presents  the  whole  list  of  comparisons  and  their 
interpretation. 
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Merlstlc  characters.     Table  2  summarizes  the  results  of  the 
analyses  of  merlstic  characters.   It  may  readily  be  seen  that,  while 
differences  related  to  the  different  environments  are  Indicated  for 
SA  and  possibly  LL,  there  are  also  differences  In  LL  In  sam.plea 
from  closely  adjacent  sink  holes  and  even  In  two  samples  from,  the 
same  sink  hole.    Thus  It  seems  that  only  SA  supports  the  separation 
of  the  two  named  subspecies. 

TABLE  2 

VARIANCE  RATIO  (F)  VALUES  FOR  COMPARISONS  OF  MERISTIC 

CHARACTERS 


Char- 

Special Series 

Main  Series 

acter 

Cl 

C8 

C, 

ClO 

m 

1.38 

0.17 

0.00 

0.51 

0.17 

0.00 

i  A 

1.87 

16. 76*** 

0.30 

0.08 

0.23 

0.93 

8.85*** 

33. 08*** 

1.04 

12.62*** 

24.91*** 

1.35 

146.94*** 

11.85*** 

15.51*** 

0.00 

1.58 

0.40 

8B 

0.06 

6.56* 

0.12 

0.04 

2.83 

0.45 

Counts  of  merlstic  characters  were  discontinued  after  analysis  of 


the  Special  Series  because  the  results  indicated  that  these  data  were  not 
of  value  in  distinguishing  the  subspecies.   During  the  course  of  my 
counting  of  scales  above  the  lateral  line,  I  observed  that  a  very  slight 
longitudinal  displacement  of  the  dorsal  fin  would  result  in  the  increase 


or  decrease  of  the  scale  count  by  one  scale.   Thus  the  variation  ob- 
served between  samples  from  the  standing -water  and  riuming- water 
habitats  might  well  have  been  due  to  variation  in  the  location  of  the 
dorsal  fin.    This  character  (UD)  does  show  significant  variation  be- 
tween the  two  habitats,  but  an  analysis  of  covariance  of  the  scale 
coimt  on  LD  did  not  reduce  the  value  of  the  variance  ratio  to  a 
nonsignificant  level  (Special  Series,  F  =  58.10***}  Main  Series, 
F  =  4.56*). 

Angle  of  gape.    Since  the  meristic  characters  were  not  to  be  in- 
cluded in  the  multivariate  analysis,  and  since  the  angle  of  gape  (AG) 
was  the  only  other  character  for  which  the  units  were  not  either 
linear  or  ratios  of  linear  units,  this  character  was  also  analyzed 
separately.    The  analysis  of  the  Special  Series  showed  that  there  was 
no  significant  variation  between  the  samples  collected  in  still  water  and 
the  one  collected  in  rimning  water  (F^  2.38)  and  that  the  variance  ratio 
for  the  comparison  of  the  two  still-water  samples  only  exceeded  the 
tabulated  F  value  for  a  5  percent  significance  level. 

For  the  Santa  Fe  Series,  one  orthogonal  comparison  (C3)  measures 
variation  between  material  collected  with  a  fifty-foot  bag  seine  and 
that  taken  with  short  seines.   Another  comparison  (C4)  compares  two 
different  collections  made  in  the  same  month,  and  a  third  (C5)  com- 
pares a  collection  made  in  February  with  one  made  in  September. 
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The  fourth  orthogonal  component  of  the  variation  (the  comparison  C^) 
was  chosen  as  the  remaining  variation  because  the  unequal  sample  sizes 
precluded  any  other  meaningful  comparison.    The  ^.-esults  of  this 
aaalyils,  which  are  presented  In  Table  18,  Indicated  no  significant 
diiferences  ia  any  of  these  components.  :    '  ^ 

In  the  Main  Series,  six  meanlngftil  comparisons  are  possible:  Cy, 
Cg,  and       are  the  same  as  they  were  for  the  meristic  characters. 
With  the  added  data  of  SGR  and  SLR,  Cxo  now  compares  the  sample 
from  Blue  Springs  with  all  three  samples  from  the  Marlon  Coxinty 
springs,  Instead  of  with  Juniper  Springs  only.    The  two  additional 
degrees  of  freedom  allow  comparisons  among  the  samples  from 
springs  In  Marlon  Co\mty:  C^  compares  the  sample  from  Jiulper 
Springs  with  the  other  two,  and  Cj2  compares  the  samples  from  Salt 
Springs  and  Silver  Glen  Springs.    In  this  analysis,  also  presented 
in  Table  18,  it  appears  that  there  may  be  some  difference  in  AG  In 
fish  from  the  two  habitats  £uid  also  some  geographic  difference.  How- 
ever, the  variation  found  in  two  samples  from  the  ssune  sink  hole  and 
the  results  of  the  analysis  of  the  Special  Series  reduce  the  confidence 
that  can  be  placed  in  this  character.   In  addition*  as  mentioned  above, 
the  precision  of  the  measurement  Involved  Is  open  to  some  question. 

Measurements .  The  comparison  Ci  for  all  measurements  in  the 
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Special  Series  indicates  that  there  are  differeaces  between  fi»h  from 
running  water  and  those  from  still  water.    For  all  bat  one  of  these 
measurements,  the  comput«d  F_  exceeds  the  0.  1  paxcent  vaiu<s  in  the 

-  table.    Comparison         which  compares  the  lime  rock  pit  sample 
with  that  from  Zamia  Sink,  indicates  significant  differences  in  eight 
of  the  twelve  measurements;  those  which  do  not  exhibit  a  significant 

value  are       HL,  SI*,  and  LAj, 
The  complete  array  of  comparisons  in  the  K^ain  Series  (Table  3) 
suggests  that  there  is  much  local  variation  in  these  characters.  In 
Jerome  Sink  and  Cow  Sink,  which  are  located  about  200  yards  apart, 
there  are  differences  in  every  character  except  HL/HW.   And  in  the 
two  samples  from  Jerome  Sink,  there  are  differences  in  all  except 
that  character  and  HL/SL.    On  the  other  hand,  the  comparison  of  all 
spring  samples  with  all  the  sink  hole  samples  (C-^)  shows  a  difference 
in  HL/HW,  but  none  in  HL/SL,  SL,  £,  or  PD.    Some  geograi^e 
variation  is  indicated  by  the  significant     values  in  Cjq  for  the  ^ 
folloA)(dng  measurements:  L,  HL,  SL,  E,  DB,  PD,  and  HL/HW.  For 
the  other  two  comparisons  in  this  series.  Juniper  Springs  fish  were 
compared  with  those  from  Silver  Glen  Springs  t^ad  Salt  Springs,  and 
the  latter  two  were  compared  with  each  other.    Between  Juniper  and 
the  other  two  springs,  there  were  differences  in  all  the  characters 
except  for  two  ratios,  HL/E  and  HL/HW.    The  Salt  Springs  and  Silver 


if 

Glen  Springe  saunples  achlblted  differences  In  all  characters  but  these 
two  and  the  ratio  PD/DB,  * 

From  the  results  of  this  analysis,  one  can  conclude  that,  with  the 
possible  exception  of  the  ratio  of  head-length  to  head- width,  none  of  the 
measurements  taken  will  allow  consistent  Identification  of  fish  from  the 
two  habitats.   Even  this  ratio  seems  to  exhibit  some  geographic  varia- 
tion. 

The  results  presented  above  Indicate  that  there  Is,  In  samples  of 
Hybopsis,  some  variation  that  Is  related  neither  to  geographic  factors 
nor  to  differences  In  the  habitat.    Since  large  collections  taken  from 
the  Santa  Fe  River,  over  a  period  of  time  and  by  various  methods, 
were  available  to  the  collections  of  the  Florida  State  Museum,  some  of 
these  samples  were  Included  as  the  Santa  Fe  Series  In  this  study.  In 
this  analysis  an  effort  was  made  to  relate  differences  In  these  m.easure- 
ments  to  differences  In  methods  or  times  of  collecting. 

With  the  five  samples  of  the  Santa  Fe  Series,  four  orthogonal  com.- 
ponents  of  the  variation  can  be  distinguished.    The  first  comparison 
Includes  the  samjple  (SFR-3)  of  25  Individuals  taken  with  short  seines 
because  this  series  was  studied  for  the  particular  purpose  of  examla> 
Ing  differences  caused  by  variable  sampling.   Inclusion  of  this  sample 
precludes  orthogonality  of  a  comparison  between  samples  taken  at  the 


same  time  in  different  years.    The  foiurth  comparison  (C^)  Is  thus  a 
measure  of  the  variation  remalaiag  after  removal  of  the  variation  du« 

to  Cj*  C4,  and  Ci^. 

The  analysis  of  C3  showed  that  there  are  significant  differences  In 
samples  taken  with  different  equipment  in  all  of  the  characters  except 
HL/HW.   In  the  comparison  of  the  two  July  samples  (C4),  only  HL/HW, 
DB,  and  PD/DB  show  significant  differences.   Between  the  February 
and  September  samples  there  are  significant  differences  Lc.  the  four 
ratios  and  in  PD  and  E  as  well.   Even  after  the  removal  of  variation 
due  to  these  three  comparisons,  there  is  still  significant  variation  in 
all  of  the  characters  except  HL./SL.    Taking  all  measurements  of  this 
series  collectively,  U  is  evident  that  there  is  no  measurement  in- 
cluded In  which  there  Is  not  some  significant  variation  among  this  set 
of  five  samples.  ' 

On  the  basis  of  the  analysis  of  variance  of  the  three  series  of 
samples,  we  may  conclude  that  there  is  no  character  among  the 
measurements  mad(;  that  will  allow  consistent  separation  of  Hybopsls 
In  Its  two  habitats.    Before  proceeding  to  the  multivariate  analysis  of 
the  twelve  characters,  it  might  be  instructive  to  select,  for  a  graphic 
analysis  such  as  that  suggested  by  Pimentel  (1953),  that  character 
which  most  nearly  allows  diagnosis  of  the  forms  in  the  two  habitats. 


If 

Graphic  Analysis  of  Angle  of  Gape 
Althoxigh,  as  was  pointed  out  earlier,  the  precision  of  the 
xneasurexnent  of  the  angle  of  gape  is  open  to  question,  and  although  it 
varies  considerably  within  one  sink  hole,  this  nevertheless  is  the  best 
single  character  for  this  purpose.    Moreover,  It  was  the  only 
character  for  which  Bartlett's  test  indicated  that  the  variance  was 
homogeneous  in  both  the  Main  Series  and  the  Special  Series,  and 
homogeneity  of  variance  is  a  prerequisite  for  Pimentel's  method. 
Figure  2  is  a  graphic  representation  of  variation  in  AG  for  all 
samples  in  the  Main  Series  and  the  Special  Series. 

It  is  obvious  in  this  chart  that  this  character  will  not  support  re- 
cognition of  the  two  subspecies  that  have  been  proposed;  nonoverlap 
of  the  entire  open  rectangle  is  required  to  indicate  subspeciflc 
differentiation  under  the  interpretation  of  the  75  percent  rule  suggest- 
ed by  Pimentel  (84  percent  of  one  popxilation  different  from  84 
percent  of  another). 
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MULTIVARIATE  ANALYSIS 
Mahalanobis'  Generali»ed  Distance  Function 

The  multivariate  technique  chosen  for  this  study  is  Mahalanobis' 
Generalised  Distance  Function  (Rao,  1952).   This  method  was  chosen 
over  Fisher's  Discriminant  Function,  Hotelling's  T^  ,  and  the  Factor 
Analysis  of  psychometry  largely  because  of  the  ease  with  which  one 
can  arrive  at  an  intuitive  interpretation  of  the  Distance  Function. 

The  Distance  Function  has  been  developed  almost  exclusively  by 
members  of  the  Indian  Statistical  Institute  of  Calcutta.   It  has  been 
used  for  distinguishing  between  neurotic  groups  (Rao  and  Slater,  1949) 
and  between  phases  of  locusts  (Nalr,  1951j  Blacklith,  1957),  as  well 
as  in  several  anthropological  investigations  (Mahalanobis,  Majumdar, 
and  Rao,  1949;  Majumdar  and  Rao,  1958).    Bailey  (1956)  examined 
the  now  classical  Partula  variation  (Cram.pton,  1916,  1932)  by  apply- 
ing this  statistic.   Although  he  questioned  the  effectiveness  of  genetic 
drift  as  a  cause.  Bailey  confirmed  Crampton's  conclusions  as  to  the 
distinctness  of  local  populations  of  this  snail. 

For  exploratory  investigations,  t^e  computation  of  the  Generalised 
Distance  (Df)has  many  advantages  over  other  multivariate  techniques. 
The  very  fact  that  Mahalanobis  chose  the  term  "distance"  for  the 
statistic  shovdd  make  it  attractive  to  students  of  speciatlon  and  racla- 
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tlon.    Some  kind  of  distance  is  certainly  involved  in  the  formation  of 
species  and  races,  although  this  distance,  regardless  of  how  it  is 
measured,  is  of  course  merely  a  measurement  of  the  probability  of 
gene  flow.    That  organisms  do  not  recognize  the  geographic  or  other 
distances  that  men  measure  is  emphasized  by  Dobzhansky's  (1951) 
definition  of  the  terms  "allopatric'*  and  "sympatrlc." 

When  a  biologist  performs  a  t_test,  he  is  in  a  sense  determining 
the  distance  between  estimates  of  the  means  of  two  populations  of 
single  variables;  he  has  converted  differences  in  morphology  to  a 
linear  scale  with  standard  units,  and  he  asks  the  question,  "Are  these 
estimates  of  the  average  values  for  this  character  far  enough  apart 
on  this  scale  to  reject  the  hypothesis  that  they  estimate  the  same 
population  mean?"     Hubbs'  and  Hubbs'  (1953)  and  Pimentel's 
(1958)  graphic  representations  of  variability  help  one  to  vlstiallze  a 
statistical  distance  between  populations  with  respect  to  a  single 
character.    With  two  uncorr elated  characters,  and  with  enough 
different  shading  techniques,  this  device  could  be  extended  by  super- 
imposing the  two  figures  at  right  angles  to  each  other.  Discordant 
variation  would,  of  course,  render  the  resulting  figure  less 
attractive  if  not  meaningless.   Further  extension  to  three  characters, 
while  difficult  to  represent  visually.  Is  theoretically  possible  by 
adding  a  third  dimension. 
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The  Generalized  Distance  Ftinction  also  assists  in  the  visualization 
of  differences  between  populations  up  to  three  dimensions.  Beyond 
that,  however,  graphic  representation  becomes  Impossible  and  the 
computations  become  laborious  even  on  an  electric  desk  calcvilator. 
However,  with  the  rapidly  Increasing  availability  of  digital  computers 
{such  as  the  IBM  650),  the  complexity  of  the  computations  should 
constitute  no  barrier,  even  to  biologists  without  mathematical 
inclination. 

Geometric  Derivation  of  D** 
In  multivariate  analysis,  geometric  Interpretation  (even  for  more 
than  three  dimeasions)  is  frequently  invoked  to  make  clearer  the 
interpretation  of  results.   For  example,  the  Generalized  Distance  be- 
tween two  populations  based  on  two  uacorrelated  characters  (Xi  and 
X^)  can  be  represented  by  the  linear  distance  between  two  points 
plotted  on  regular  Cartesian  coordinates.    The  formula  for  this 
Distance  is: 


(X,i  -  Xu)^  -  Xz;)^ 


where  Xji  and  Xiz  are  the  means  of  XI  for  the  first  and  second  popu- 
lations, respectively;  X21  and  X22  ttie  means  of  X2;  and  sl^  and 
the  pooled  estizxuites  of  the  variances  of  Xx  and  X2.   Thus,  the  use  of 
the  Generalized  Distance  requires  the  assiimptlon  that  the  variance  of 
each  variable  is  the  same  in  all  populations.    This  assumption  may 


not  always  be  valid,  but  Rao  (1952)  has  stated  that  departures  from  it 
are  not  serious  if  the  measures  are  at  all  related  and  if  s^  is  es- 
timated from  a  large  number  of  measurements.   For  operations  with 
this  statistic,  extraction  of  the  square  root  to  obtain  D  Is  generally 
delayed  until  the  last  operation,  X>  being  used  directly.   This  prac- 
tice will  be  followed  in  the  present  analysis. 

If  this  Xi  and       are  not  uncorrelated,  Cartesian  coordinates  are 
not  valid  for  representation  of  the  distance.    In  this  situation,  coordi- 
nates are  constructed  which  are  not  perpendicular  but  which  are  in- 
clined to  each  other  at  an  angle  whose  cosine  is  r,  the  coefficient  of 
correlation  between  tiie  variables.   In  the  special  case  of  uncorrelated 
variables,     =  0,  and  the  angle  whose  cosine  is  sero  is  90^,  so  that 
the  axes  are  perpendlciilar.    The  formula  given  above  must  be 
modified  to  include  any  lack  of  orthogonality  due  to  correlated 

variables.   In  algebraic  notation:  ^ 
(X„  ~%<zf  ^   (Xh  -Xi2)(Xii  ->?2Z)  ^(X^l 

D  =  — +  t:,  ^ 

with  the  symbols  as  above  except  that  B12,  is  the  co variance  of  var- 
iables Xi  and  X^. 

It  is  evident  that  such  formulae  will  become  increasingly  unwieldy 
as  the  number  of  variables  increases,  so  It  Is  much  more  convenient 
to  use  summation  rather  than  algebraic  notation.  Letting 


di=  (Xii-Xi2).  d2=  (X2i-X22)--l.e.,  djc  (Xq-Xu)  and  sij  be  the 
covariance  of  the       and  j*^  variables  (or  the  variance  of  the  i^ 
variable  when  I  =  j): 

The  superBcript  indicates  the  inverse  of  the  variance -covariance 

matrix.    Both  Rao  (1952)  and  Roycc  (1957)  present  methods  for 
simplifying  the  computations  of  D  . 

Thus,  we  can  see  that  by  scaling  Cartesian  coordinates  in  stan- 
dard deviation  units  and  applying  the  Pythagorean  theorem,  the 
distance  between  population  means  of  the  two  independent  variables 
can  be  estimated;  and  that  this  concept  can  be  extended  to  correlated 
variables  and  to  three  or  even  more  dimensions.   Although  the 
computations  would  be  involved  and  laborious ,  the  use  of  digital 
computers  and  the  ease  of  Interpretation  make  Mahalanobls' 
Generalized  Distance  Function  an  extremely  useful  statistic  for 
simultaneous  suialysls  of  several  taxonomlc  characters.   It  is  also 
Interesting  to  note  at  this  point  that,  contrary  to  Mayr,  Llnsley.  and 
Uslnger's  (1953)  assertion  that  the  coefficient  of  difference  had  not 
been  extended  to  multiple  characters,  D    Is  one  of  several 
statistically  valid  measures  of  differences  In  multiple  characters 
which  antedate  their  D.D.  by  about  a  quarter  of  a  century. 


34 


Computation  of  D 
Since  there  is  not  yet  available  a  program  for  computing  on 
the  IBM  650,  this  computation  was  performed  using  parts  of  several 
other  programs.    In  order  to  compute  also  the  single  factor  statistics 
here  reported,  the  sum  of  squares  and  products  matrix  was  determin- 
ed separately  for  each  sample.    Then,  because  the  variances  were 
not  the  same  in  all  samples  and  because  Bao  (1952)  suggested  that 
this  would  not  invalidate  the  technique  if  the  variances  and  co- 
variances  were  based  on  a  large  number  of  specimens,  the  sums  of 
sq\iares  and  products  for  all  of  the  635  specimens  were  pooled  for 
an  estimate  of  the  variance-covariance  matrix  based  on  620  degrees 
of  freedom.    This  matrix  was  then  inverted;  the  differences  between 
the  means  of  the  12  measurements  for  all  45  pairs  of  samples  in  the 
Main  Series  and  Special  Series  were  determined,  and  the  vector- 
matrix-vector  miiltiplication  was  performed  to  arrive  at  a  value  of 

2  y 
D   for  each  pair  of  samples.    Table  4  presents  these  values  of  D  . 

The  variance-covariance  matrix  and  its  inverse  are  presented  in  the 

Appendix  in  Tables  10  and  11,  respectively. 

Interpretation  of 

If  the  two  subspecies  described  by  Hubbs  and  Crowe  (1956)  are 

valid,  a»d  if  this  divergence  is  reflected  in  the  morphometric 
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characters  included  in  tluB  analysis,  the  populations  from  which  these 
samples  were  taken  would  form  two  distinct  clusters.    V/hilc  the 
criteria  for  clusters  are  not  rigidly  fixed,  "any  two  groups  belonging 

to  the  same  cluster  should,  at  least  on  the  average,  show  a  smaller 

2 

D   than  those  belonging  to  different  clusters"  (Hao,  1952). 

■    ■    TABLE  4 
VALUES  OF  D^.  CORRECTED  FOK  BIAS 


ZMS 

SFR 
-1 

HPS 

BLR 

SGR 

SLR 

JUR 

CWS 

JRS 
-1 

SFR-1 

6. 1 

HPS 

2.8 

2.0 

BLR 

5.5 

3.7 

3.0 

SGR 

5.5 

2.4 

2.6 

33.6 

SLR 

33.4 

6.4 

4.5 

37.4 

3.4 

JUR 

7.0 

2.2 

5.8 

162.9 

4.  1 

7.0 

CWS 

4.0 

3.8 

1.3 

157.7 

5.0 

8.3 

9.1 

JRS-1 

7.5 

7.5 

7.2 

95.9 

20.3 

47.5 

16.8 

16.8 

JRS-2 

8.8 

9.3 

6.9 

54.4 

33.6 

6.4 

398.2 

392.3 

7.4 

The  first  question  in  the  present  case  is,  therefore,  whether  two 
clusters  can  be  found,  one  consisting  of  the  running-water  samples 
and  the  other  comprising  the  still-water  samples.    When  such 
clusters  are  formed,  the  intracluster  average       is  26.  31  for  the 
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running  water  and  45. 17  for  the  standing  water.    The  average 
Intercluster        is  36.92;  thus,  this  arrangement  does  not  fulfill 
the  criterion  suggested  above. 

Since  the  habitats  of  the  subspecies  in  the  original  description  are 
not  delimited  as  above,  but  rather  as  "springs  and  spring-fed  creeks" 
contrasted  with  "sinks  and  caves"  (Hubbs  and  Crowe,  1956),  U 
might  be  necessary  to  leave  SFR-1  and  HPS  out  of  consideration  and 
cluster  the  other  samples  according  to  these  habitats.    With  this 
arrangement,  the  spring  intracluster  average        is  41.40;  that  of 
the  sink  hole  habitat  is  72.27;  and  the  average  d£_  between  the  two 
habitats  is  54.82.    This  arrangement  likewise  falls  to  satisfy  the 
criterion  for  clustering. 

Since  neither  of  these  a  priori  cluster  formations  will  satisfy  the 
criterion,  an  alternative  procedure  Is  to  search  for  clusters  In  the 
manner  suggested  by  K.  D.  Tocher  to  Rao  (1952).    Beginning  with  two 
samples  with  a  low       ,  other  samples  are  added  to  the  cluster  imtil 
such  aua  addition  causes  a  large  Increase  either  in  the  average  or 
in  the  average  increase  in  Lr.     The  simultaneous  development  of 
two  clusters  sometimes  allows  a  better  decision  regarding  the 
retention  or  discarding  of  a  sample  in  a  cluster. 

Following  this  procedure,  and  beginning  the  two  clusters  with  ZMS 
and  SFR-1  of  the  Special  Series,  the  following  clusters  are 
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dlBcernible: 

Ai,  including  samples  SFR-1,  JUR,  SGR,  and  SLR; 

Az»  including  samples  ZMS  and  CWS; 

A^t  includiiig  samples  JRS-1  and  JRS-2;  and 

A4,  including  only  sample  BLR. 
Sample  HPS  is  not  incliaded  in  any  of  these  clusters,  because  the 
value  of       between  this  and  each  of  the  other  samples  is  so  low 
that  HPS  can  be  added  to  any  of  the  clusters  at  any  time  without 
violating  the  stipulation  for  cluster  formation.    The  average 
Intercluster  and  Intracluster  values  for       are  presented  in  Table  5. 

TABLE  5 

INTRACLUSTER  AND  INTERCLUSTER  VALUES  FOR  Df  IN  THE 
CLUSTERING  CONFIGURATION  "A" 


A| 

A2 

A$ 

M 

HPS 

Al 

4.25 

9.78 

64.38 

59.4 

4.05 

A2 

9.78 

4.0 

132.82 

81.6 

4.70 

A3 

64.38 

132.82 

7.4 

75.15 

7.17 

A4 

59.4 

81.6 

75.15 

3.0 

HPS 

4.05 

4.70 

7. 17 

3.0 

Cluster  Al  includes  samples  SFR-1,  JUR,  SGR,  and  SLR 
Cluster  Az  includes  samples  ZMS  and  CWS 
Cluster  A3  includes  samples  jRS-1  and  JRS-2 
Cluster  A4  includes  only  sample  BLR 
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Obviously  this  configuration  does  satisfy  the  criterion  of  lower 
average  values  for  the  Intracluster       ;  It  Is  hereafter  designated  as 
"A"  and  Is  Illustrated  in  Figure  3.   Clusters  Aj  and       seem  to  be 
much  closer  together  thzm  the  others,  and  this  stiggests  that  at  least 
one  of  the  members  of  A2  might  belong  In  Aj.   Addition  of  CWS  to 
Aj  Increases  the  average  Intracluster       to  5.17,  which  Is  still 
considerably  below  the  Intercluster  values  of  D*^. 
•   Thus,  an  alternative  configuration  may  Inclxxde  CWS  with  Aj,  and 
•ince  the       between  BIJl  and  ZMS  Is  less  than  that  between  JRS-l 
and  JRS-2  and  Is  very  little  more  than  the  Intracluster  average 
In  Ap  these  two  samples  are  also  combined  to  make  a  cluster.  This 
configuration  Is  designated  "B"  and  Is  also  Illustrated  In  Figure  3, 
The  Intracluster  and  Intercluster  values  for       are  presented  m 
Table  6. 

While  It  was  difficult.  If  not  Impossible,  to  discern  any  pattern 
In  the  variation  of  these  twelve  characters  by  the  methods  of  single 
factor  analysis,  some  sort  of  structure  seems  to  emerge  In  either 
one  of  the  cluster  configurations  presented  above.   Considering  only 
tlie  criterion  of  maximum  Intercluster  and  minimum  Intracluster 
values  for  D^,  the  arrangement  In  "B"  seems  preferable,  but 
configuration  "A"  also  has  some  points  to  recozxmiend  It.   In  the  first 
place,  with  "A"  most  of  the  samples  from  the  stream  habitat  form  a 
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TABIuE  6 

XNTRAC LUSTER  AND  INTERC LUSTER  VALUES  FOR  Df  IN  THE 
CLUSTERING  CONFIGURATION  "B" 


»2 

^3 

HPS 

5.17 

45.13 

91.98 

3.26 

45.13 

5.5 

16.66 

2.85 

% 

91.98  , 

16.66 

7.4 

7.05 

HPS 

3.26 

2.85 

7.05 

Cluster  Bi  includes  samples  SFR-1,  JUR,  SCR,  SLR, 
and  CWS 

Cluster  B2  includes  samples  ZMS  and  BLR 
Cluster  B3  includes  samples  JRS-1  and  JRS-2 


distinct,  separate  cluster.    The  failure  of  BLR  to  be  included  in  this 
cluster  can  be  eiqplained  because  this  sample  was  taken  in  the  back- 
water aroiind  the  boil  of  the  spring  during  a  period  of  very  high  water. 
Therefore,  it  may  have  sampled  a  population  which  is  less  similar  to 
the  spring  form,  either  because  of  actual  relationship  to  the  sink  hole 
populations  or  because  of  a  sorting  out  in  this  area  of  less  positively 
rheotropLc  individuals. 

Also,  In  "A"  It  appears  that  there  may  be  two  major  axes  In  the 
distribution  of  the  clusters:  one  placed  horizonally  in  the  figure  and 
the  other  running  almost  vertically  near  Aj  and  A2.    Under  a 
hypothesis  of  such  axes,  the  horizontal  one  might  be  Interpreted  as 
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representing  variation  related  to  the  sink  hole  habitat  and  the  vertical 
one  as  related  to  the  stream  habitat.    This  would  thus  lead  to  the 
conclusion  that  the  range  of  variability  of  the  sink  hole  form  is  so 
broad  that  it  extends  to  both  sides  of  the  stream-form  cluster.  This 
could  also  account  for  the  small  Distances  which  allow  inclusion  of 
CWS  in  Cluster  Aj  and  for  the  fact  that  HPS  can  be  added  to  any  of 
the  clusters  without  increasing  the  intracluster  Distance  excessively. 

Although  the  configuration  in  "B"  meets  the  criterion  for  cluster 
formation,  it  is  possible  that  this  figure  is  too  greatly  oversimplified. 
This  arrangement  does  emphasize  the  possibility,  mentioned  above, 
that  BLR  Is  a  sample  more  representative  of  the  sink  hole  form  than 
of  the  spring  form.    This  hypothesis  is  further  supported  by  the  fact 
that  Cluster  B2  is  closer  to  B3  than  it  is  to  Bj. 

There  is  actually  little  or  no  contradiction  in  the  two  figures.  If 
"A"  were  to  be  rotated  until  BLR  were  projected  onto  A2,  a  figure 
very  much  like  "B"  would  be  obtained.    The  configuration  in  "A"  is 
to  be  preferred  because  this  rotation  would  in  effect  bring  about  a 
greater  reduction  in  the  information  conveyed  by  the  figure.    The  en- 
tire process  of  cluster  formation  that  produced  "A"  resulted  in  a  loss 
of  information,  but  these  rotations  were  necessary  to  obtain  a  struc- 
ture that  could  be  interpreted  and  graphically  represented. 


DISCUSSION  AND  CONCLUSIONS 
In  recent  years  the  nonrandomness  of  samples  has  come  gradually 
to  be  seen  as  a  major  source  of  error  in  taxonomic  work.  Royce 
(1957)  commented  on  the  difficxilty  of  insuring  randomness  in  samples 
of  fish,  and  Plmentel  (1958)  urged  systematlsts  to  seek  complete 
raadomia*tlon  of  collections.  A  statistical  test  in  which  departures 
from  randominess  are  not  serious  Is  much  to  be  desired  because  most 
taxonomic  work  and  studies  of  variation  are  forced  to  rely  to  a  large 
extent  on  museum  collections.    Plmentel  suggested  the  application  of 
a  test  for  homogeneity  of  variances,  this  to  be  followed  by  transfor- 
mations of  the  variables  until  homogeneity  was  obtained.  An 
alternative  procedure  Is  the  use  of  Mahalanobls'  D^  ,  which  can  be 
corrected  for  bias  due  to  different  sample  nximbers  and  which.  If  the 
estimates  of  the  variances  are  based  on  a  large  number  of  degrees  of 
freedom,  Is  not  greatly  affected  by  departures  from  homogeneity  of 
variances. 

The  material  on  which  this  study  Is  based  Is  subject  to  all  the 
departures  from  randomness  to  be  fovind  in  museum  collections. 
Random  samples  from  such  widely  divergent  habitats  as  caves,  sink 
holes,  and  springs  would  be  practically  impossible  to  obtain,  since  the 
same  collecting  techniques  are  not  equally  effective  In  these  sltua- 
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tions.   For  example,  on  a  day  when  fifteen  minutes'  trapping  effort 
yielded  more  than  30  fish  in  Jerome  Sink,  an  hour's  trapping  in  Blue 
Springs  yielded  none,  even  though  great  niunbers  of  fish  could  be 
observed  swimming  about  the  trap  and  bait. 

Thus,  even  though  the  results  must  be  Interpreted  with  caution  be- 
cause of  the  heterogeneity  of  variance,  the  analysis  of  variance  tech- 
nique failed  to  reveal  any  characters  for  which  variation  could  be  re- 
lated to  habitat  distribution.    Furthermore,  analysis  of  a  series  of 
samples  from  one  area  of  the  Santa  Fe  River  showed  that  most  of 
these  variables  are  subject  to  significant  sampling  variation. 

However,  when  the  data  were  processed  through  a  multivariate 
technique  with  D^,  a  structure  emerged  which  is  generally  in  accord 
with  the  results  hypothesised  above  from  the  distribution  of  Hybopsis. 
There  are,  of  course,  some  rather  large  departures  from  the 
predicted  results,  but  these  can  be  explained  on  the  basis  of 
deviations  either  in  the  population  structui;e  or  in  the  sampling 
procedure.   Some  of  these  deviations  might  also  have  affected  the 
single  factor  analysis,  but  in  the  absence  of  a  basic  pattern,  the 
deviations  could  not  be  identified. 

The  reduction,  by  clustering,  of  the  twelve-dimensional  hyper- 
space  in  which  the  samples  were  distributed  by  the  Generalized 
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Distance  analysis  resiilts  In  a  distribution  of  the  clusters  along  two 
znajor  axes  (Figure  3).    The  horizontal  axis  in  this  figure  can  be 
taken  as  representing  the  range  of  variation  of  the  sink  hole  form  of 
Hybopsls  harperl.     Hypothetlcally,  a  single  Isolate  from  a  sink  hole 
popxilatlon  would  produce  a  subpopulatlon  from  which  a  sample  could 
be  found  somewhere  along  tills  axis.   A  second,  less  well-defined 
axis  can  be  distinguished,  and  this  one  may  represent  the  dimension 
through  which  the  stream  form  varies.   On  this  axis,  BLR,  a  sample 
from  a  part  of  a  spring  with  much  reduced  current,  is  found  near 
one  extreme.    The  other  stream  samples  form  a  distinct  cluster 
close  to  the  intersection  of  this  axis  and  the  sink  hole  axis. 

To  continue  this  geometric  Interpretation:  we  can  visualize  a 
change  In  morphometrlc  characters,  regardless  of  cause,  to  be 
manifested  as  movement  of  a  population  through  a  hyperspace  which 
Is  represented  here  by  only  two  or  tiiree  dimensions.  The 
characteristics  of  the  stream  habitat  produce  movement,  however 
slight,  mainly  In  the  direction  of  the  A1-A4  axis,  and  the 
characteristics  of  the  sink  hole  habitat  produce  changes  along  the 
A2-A3  axis.    Since  sink  hole  populations  are  almost  certainly 
derived  ultimately  from  stream  populations,  the  sink  hole  samples 
may  be  considered  to  have  moved  away  from  the  cluster  of  stream 
popiilatlons  along  the  axis  of  the  sink  hole  habitat. 
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The  extreme  position  oi  BLR  on  the  stream  habitat  axis  has 
already  been  mentioned  and  a  possible  explanation  suggested.  There 
are  several  other  points  which  require  further  explanation:  the 
location  ol  CWS  which  would  allow  its  inclusion  in  the  stream  sample 
cluster  (A^)  and  the  very  large  distance  between  tide  sample  and 
JRS-2.  while  ita  Distance  from  JRS-1  is  much  less;  the  small 
Distances  between  HPS  and  all  the  other  samples;  and  the  relatively 
large  Distance  between  JRS-1  and  JRS-2.  ' 

From  1955  until  1959.  the  general  water  table  was  low,  and  in  Cow 
Sink  only  about  fifteen  sqtiare  feet  of  water  surface  was  exposed. 
During  this  time  the  population  of  Hybopsis  was  very  small.  In 
fact,  on  several  occasions  one  could  see  only  two  or  three  fish  in  the 
sink  hole.   Throughout  this  time.  Cow  Sink  was  used  as  a  source  of 
irrigation  water  during  the  tobacco  growing  season.   In  1959  tine 
water  level  rose  considerably,  and  the  Hybopsis  population  of  Cow 
Sink  increased  greatly.    Two  factors  in  this  situation  might  explain 
-  the  small  distance  between  CWS  and  the  cluster  of  stream  samples. 
First  the  environment  is,  with  respect  to  the  fundamental  factor  of 
water  flow,  comparable  to  a  small  intermittent  spring  with  a  flow  of 
about  10,  000  gallons  per  hour  (the  output  during  the  irrigation  season). 
Thus,  some  enviroiunentally  produced  phenotypic  eftect  might  hold 
CWS  nearer  the  intersection  of  the  two  main  axes.    On  the  other 


]iand»  the  population  which  is  ssusipled  by  CWS  might  be  derived 
either  from  the  very  reduced  population  that  survived  the  period  o£ 
low  water  or  from  a  colonization  by  new  migrants  into  the  sink  hole. 
In  either  case,  the  genetic  composition  of  the  newly  increased 
population  would  be  similar  to  that  of  the  stream  population.  Ob- 
viously, new  migrants  would  be  expected  to  bring  a  random  sample  of 
the  genes  of  the  population  from  which  they  came.   And  if ,  as 

-'    V  -a' 

Lerner  (1954)  has  suggested,  "natural  selection  will  tend  to  favor 
organisms  clustering  arovmd  mean  values  for  all  characters,"  the 
survivors  of  the  period  of  reduced  population  would  be  expected  to 
have  a  genetic  constitution  like  that  of  the  average  population,  which, 
according  to  Figure  3,  is  the  stream  form.   In  any  case,  even  though 
CWS  Is  Included  in  the  cluster  with  ZMS,  it  is  much  closer  to  the 
stream -form  cluster  than  is  ZMS. 

The  peculiar  relationship  of  the  Distances  between  JRS-1,  JRS-2, 
and  CWS  is  also  related  to  the  physical  changes  described  above,  as 
is  the  relatively  large  Distance  between  jRS-1  and  JRS-2.  Be- 
ginning around  1956,  Jerome  Sink  also  was  used  for  irrigation. 
Therefore,  the  sample  taken  in  1959  (JRS-1)  was  from  a  population 
subject  to  the  same  conditions  as  those  described  above  for  CWS, 
except  that  at  no  time  was  the  population  of  Jerome  Sink  so 
drastically  reduced.    This  population  may  have  thus  moved  along 
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tiie  sink  hole  axis  back  toward  the  stream-form  cluster  after 
irrigation  began.   By  the  same  token.  Irrigation  between  the  times 
when  JRS-2  and  JRS-1  were  collected  would  have  reduced  the 
effectiveness  of  Isolation  and  have  allowed  the  Distance  between  the 
population  samples  by  CWS  and  that  sampled  by  JRS-1  to  decrease. 

The  situation  of  HPS  is  more  difficult  to  explain.    The  greatest 
Distance  between  HPS  and  any  other  sample  is  that  between  HPS 
and  JRS-1  {0^=7.2);  this  value  Is  less.than  that  between  JRS-1  and 
JfiS-2.   The  population  sampled  by  HPS  definitely  represents  a  new 
Isolate  because  the  pond  was  formed  by  scooping  out  a  mud  hole  In 
the  bottom  of  the  lime  rock  mine.   Either  by  influx  of  many  migrants, 
or  by  the  very  rapid  reproduction  of  a  few,  the  population  increased 
to  an  estimated  several  thousand  In  a  period  of  two  years.   Again  It 
Is  possible  to  Invoke  the  idea  of  natural  selection  favoring  the  • 
average  to  account  for  the  small  distance  between  HPS  and  the 
stream-form  cluster,  but  neither  this  mechanism  nor  the  hypothesis 
of  large  scale  tmlgratlon  will  explain  the  sanall  distance  between 
HPS  and  the  sink  hole  samples.   In  fact,  the  only  hypothesis  that 
seems  to  explain  the  location  of  HPS  Is  that  the  movement  through 
the  hyper  space  was  not  limited  to  either  of  the  two  major  axes  herein 
described.    Nor  was  it  limited  to  movement  along  both  these  axes. 
The  movement  of  this  population  was  apparently  in  some  direction 
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such  that  projection  to  these  axes  places  HPS  near  each  of  the 
clusters.   This  condition  could  be  realised  If  the  aspects  of  this 
habitat  that  are  similar  to  the  stream  habitat  (e.g. ,  direct  sunlight, 
high  basic  productivity)  produce  movement  that  can  be  referred  to 
that  axis ,  and  if  those  aspects  similar  to  the  sink  hole  habitat 
(e.g. ,  no  current,  no  submerged  or  emergent  tracheophytes)  pro- 
duce movement  that  can  be  referred  to  that  aucls. 

Thus,  the  use  of  Mahalanobls'        and  the  clustering  technique 
allow  formation  of  an  array  of  samples  from  the  various  habitats 
that  is  meaningful  because  it  raises  fundamental  questions  concern- 
ing the  genetic  and  ecological  relationships  among  the  popxilations. 
Moreover,  this  Interpretation  of       provides  a  model  lor  the 
morphometric  divergence  of  subpopulations  which  has  previously 
been  Interpreted  as  subspeciatlon. 

A  final  conclusion  to  be  drawn  from  this  analysis  is  that  the 
structure  of  the  clusters  of  samples  does  not  warrant  recognition 
of  two  subspecies  in  any  conventional  taxonomic  sense.  The 
populations  inhabiting  springs  and  spring-fed  streams  are  similar 
enough  to  each  other  for  consideration  as  a  single,  discrete  tax- 
onomic entity.    On  the  other  hand,  the  divergence  of  the  sink  hole 
populations  is  not  consistently  In  one  direction,  as  one  would  expect 


It  to  be  If  there  were  two  atnd  oixly  two  subspecies  of  Hybopsis 
harjperi.   The  alternatives  that  present  themselves  are:  (1)  to  name 
and  describe  several  subspecies  from  sink  holes,  producing  even 
more  discontinuities  aiad  interdlgitations  in  aooaps  of  the  ranges  of 
the  subspecies,  or  (2)  simply  to  describe  H.  harper  I  as  a  variable 
species,  In  which  local  Isolates  diverge  from  the  typical  stream 
form  In  at  least  two  directions.    I  prefer  the  latter  course  of  action. 


SUMMARY 

Single  factor  analysis  for  five  meristic  characters  of  450  fish  in 
eight  samples  of  Hybopsia  harperi  reveal  little  variation  in  these 
characters  in  the  various  subpopulations .    The  character  most 
strongly  supporting  the  validity  of  the  subspecies  now  recognized 
taxonomically  is  the  nimiber  of  scales  above  the  lateral  line. 

Similar  analysis  of  eight  linear  measurements  and  four  ratios  for 
those  450  fish  and  for  90  more  in  three  other  samples  failed  to  show 
any  pattern  of  variation  that  could  be  related  to  the  habitats  Varia- 
tion in  the  angle  of  gape  was  also  studied;  and  although  there  seemed 
to  be  a  consistent  difference  between  samples  from  the  two  different 
kinds  of  habitat,  the  results  of  graphic  analysis  of  this  character  and 
the  questionable  precision  of  the  technique  of  measurement  preclude 
its  use  as  a  diagnostic  subspecific  character  in  H.  harperi. 

To  investigate  sampling  variation,  the  measurements,  ratios, 
and  angle  were  analyzed  in  a  series  of  samples  from  the  same 
location.    This  analysis  indicated  ttiat,  with  the  possible  exception  - 
of  the  angle,  any  one  of  these  characters  is  likely  to  show:  (1)  dif« 
ferences  in  samples  collected  by  different  methods;  and  (2)  dif* 
ferences  due  to  unidentified  causes* 

After  it  was  demonstrated  that  single  factor  analysis  of  the 
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material  available  for  this  study  did  not  reveal  any  characters  with 
concordant  variation*  the  measurements  and  ratios  were  subjected 
to  multivariate  analysis  with  MahalanobU'      .     Through  use  of 
this  technique  variational  structure  does  emerge »  atnd  the  estimated 
Distances  between  the  various  populations  permit  the  following  - 
conclusions  concerning  the  relationship  between  the  poptilatlons: 

1.  In  the  twelve -dimensional  hyperspace  formed  by  this  analysis, 
there  Is  a  major  axis,  along  which  the  sink  hole  populations  vary 
from  each  other. 

2.  Intersecting  this  axis  there  may  be  another  which  represents 
the  range  of  variation  of  the  spring  and  stream  populations. 

3.  At  or  near  the  point  of  intersection  of  these  axes,  the  stream 
populations  form  a  distinct  cluster. 

4.  A  model  for  morphometrlc  chamge  Is  suggested.  In  which 
change  is  Interpreted  as  movement  along  the  axis  corresponding  to 
the  habitat  of  an  isolated  population. 

5.  Departures  from  the  above  pattern  can  be  explained  on  the 
basis  of: 

(a)  changes  In  the  physical  environment, 

fb)  rapid  Increase  In  popiUatlon  numbers  which  allows 

survival  of  extreme  variants,  or 
(c)  selection  for  the  average  during  periods  of  lowered 
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water  level  and  reduced  population  (genetic 
homeostasis).  ... 
6.    The  failure  of  the  sink  hole  populations  to  form  one  distinct 
cluster  and  the  fact  that  movement  along  the  "sink  hole  axis"  may  be 
in  either  direction  seem  to  be  grounds  for  concluding  that  the 
taxonomic  separation  of  Hybopsis  harperi   into  twot  ajid  only  two, 
subspecies  does  not  satisfactorily  reflect  the  natural  situation. 
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This  appendix  contain*  the  complete  statistical  analyses  referred  to 
in  the  text.   In  most  instances,  the  entries  in  the  tables  contain  fewer 
digits  than  were  actually  used  in  the  computations.    For  calculating 
means,  sums  of  squares  and  products,  and  standard  deviations,  ten 
digits  were  used.    For  Table  7,  the  means  and  standard  deviations  of 
the  measurements  were  converted  to  millimeters. 

The  matrix  operations  involved  in  computing  D    were  performed  in 
floating  point  arithmetic;  therefore,  the  numbers  actually  used  con- 
tained eight  digits  in  scientific  notation.   For  example: 

31,768.453  =  3.  1768453  x  10*. 
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TABLE  7 

MEANS  AND  STANDARD  DEVIATIONS 
FOR  ALL  CHAKACTEKS  IN  ALL  SAMPLES 


AG 

L 

DB 

HL 

SL 

m ' 

LD 

ZMS 

27. 8 

31.6 

6.  9 

7.  9 

2.  1 

2.  3 

16.  2 

\4.  00) 

IL  AB.\ 

(o.  45) 

(1.  52) 

(1.43) 

(0.  50) 

(0.  31) 

(3.  33) 

HPS 

26.5 

31.9 

8.  1 

8.0 

2.2 

2.2 

16.  1 

(3.  74) 

(4. 19) 

(1.02) 

(1.03) 

(0.  34) 

(0.23) 

(2.  26) 

SFR-1 

27.9 

37.6 

8.0 

8.6 

2.4 

2.7 

19.2 

(3.46) 

(4.82) 

(1.02) 

(1.01) 

(0.  34) 

(0.  24) 

(2.53) 

SFR-2 

26.6 

37.5 

8.8 

8.2 

2.2 

2.5 

18.7 

(3.25) 

(6.20) 

(1.33) 

(1.34) 

(0.41) 

(0.  30) 

(3.  32) 

SFR-3 

29.7 

33.6 

8.0 

7.  1 

1.9 

2.2 

17.0 

(4.00) 

(3.  33) 

(0.97) 

(0.63) 

(0.25) 

(0.  16) 

(1.68) 

SFR-4 

28.9 

32.9 

7.3 

7.0 

1.8 

2.2 

16.0 

(3.25) 

(2.59) 

(0.86) 

(0.63) 

(0.25) 

(0.  15) 

(1.41) 

srR-5 

28.  1 

31.3 

6.6 

7.0 

1.8 

2.3 

15.3 

(2.69 

(3.  19) 

(0.72) 

(0.79) 

(0.23) 

(0.22) 

(1.65) 

SFR-6 

27.5 

38.  1 

8.6 

9.6 

2.2 

2.4 

18.7 

(2.  36) 

(3.09) 

(0.91) 

(0.76) 

(0.  16) 

(0.27) 

(1.86) 

JRS-1 

31.8 

34.6 

8.5 

7.7 

1.9 

2.3 

18.5 

(5.47) 

(7.44) 

(1.89) 

(1.53) 

(0.48) 

(0.  34) 

(3.80) 

JRS-2 

28.4 

23.6 

5.  3 

5.5 

1.3 

1.7 

12.5 

(3.46) 

(3.42) 

(0.79) 

(0.75) 

(0.25) 

(0.  16) 

(1.62) 

CWS 

28.5 

36.3 

8.0 

8.7 

2.3 

2.5 

18.8 

(5.57) 

(7.79) 

(2.31) 

(1.78) 

(0.58) 

(0.31) 

(4. 46) 

BLR 

29.7 

28.4 

6.2 

6.3 

1.6 

2.  1 

14.6 

(2.83) 

(4.86) 

(1.02) 

(0.98) 

(0.  33) 

(0.  24) 

(2.51) 

SGR 

25.3 

32.5 

7.0 

7.4 

1.8 

2.3 

15.9 

(3.00) 

(4.94) 

(1.  17) 

(0.95) 

(0.  27) 

(0.  19) 

(2. 46) 

SLR 

23.9 

24.4 

5.4 

5.7 

1.4 

1.9 

12.4 

(3.00) 

(6.02) 

(1.14) 

(0.  35) 

(0.  35) 

(0.27) 

(2.88) 

JUR 

25.9 

36.8 

6.9 

8.2 

2.3 

2.6 

18.5 

(3.60) 

(5.74) 

(1.  17) 

(1.17) 

(0.41) 

(0.29) 

(2. 84) 
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TABLE  7- -Continued 


PD 

HW 

HL/SL 

T  T  T     /  TT* 

HL/£ 

HL/HW 

A  2 

•P.  09 

3-  35 

1.  91 

0.  37 

\<J»  1  91 

(0.  41) 

(0.  28) 

(0.  12) 
I"*  •••1 

(0.  04) 

UT3C 

A  A 

7< 

•7*  3^ 

1.81 

;  1 0.  33 

/O*  25) 

(0.  19) 
\w.  *  7; 

(0.  12) 

(0.03) 

or  K-  1 

A  k 

J*  WW 

1.  98 

0.  36 

to  (iO\ 

/O.  20) 

(0.  21) 

(0. 08) 

(0.03) 

A  A 

71 

3-  24 

1.85 

0.  33 

to  71\ 

/O.  27) 

to.  25) 

(0. 06) 

(0.02) 

c  irij  ^ 
ai?  IV-  5 

7  ft 

3.  16 

1.80 

0.  36 

/ft  2<J\ 

<0  43\ 

(0.  30) 

(0.  15) 

(0.08) 

(0.02) 

7  7 

%  ft 

•7.  O 

J*  7*0 

3.  18 

1. 85 

0.  36 

(0.  34) 

to.  23) 

(0.  20) 

(0.07) 

(0.02) 

J*  t 

3.  93 

3.  10 

1. 90 

0.  39 

(0.40) 

(0.  19) 

(0.  16) 

(0.06) 

(0.03) 

or  i\-o 

3.  97 

J.  7  • 

3.  55 

1.  90 

0.  37 

to.  27i 

(0.  351 

(0.  25) 

(0.25) 

(0.08) 

(0.03) 

J  Jew-  * 

4.  3 

*X*  *7 

4. 06 

3. 36 

1.80 

0.32 

(0.53) 

(0.96) 

(0.39) 

(0.  33) 

(0.  14) 

(0.03) 

JRS-2 

1.8 

3.0 

4. 15 

3.  14 

1.83 

0.34 

(0.29) 

(0.41) 

(0.  34) 

(0.  24) 

(0.08) 

(0.02) 

cws 

2.8 

4.8 

3.80 

3.49 

1.82 

0.36 

(0.62) 

(1.07) 

(0.  38) 

(0.  39) 

(0.  13) 

(0. 04) 

BLR 

2.2 

3.5 

4.01 

3.07 

1.83 

0.36 

(0.  36) 

(0.58) 

(0.  36) 

(0.  17) 

(0.07) 

(0.02) 

SGR 

2.5 

3.9 

4.05 

3.  18 

1.91 

0.36 

(0.  35) 

(0.56) 

(0.  19) 

(0.20) 

(0.09) 

(0.02) 

SLR 

2.0 

3.0 

4.  24 

3.06 

1.95 

0.  36 

(0.46) 

(0.61) 

(0.  30) 

(0.  20) 

(0.07) 

(0.02) 

JUR 

2.9 

4.3 

3.59 

3.  18 

1.91 

0.34 

(0.41) 

(0.57) 

(0.25) 

(0.  25) 

(0.  12) 

(0. 02) 

'  * 

.   .             ■  i 

\    »  ■   ■  ', 

•'■  '* 

•.  ••>•  '  .' 

TABLE  7- -Continued 


D 

A 

C  A 

ZMS 

8.  00 

7.88 

35.8 

5.21 

2.99 

(0.  14) 

(0.  38) 

(1.48) 

(0.43) 

(0.  27) 

HPS 

7.99 

8.07 

36.7 

5.41 

3.11 

(0.22) 

(0.  38) 

(1.01) 

(0.49) 

(0.  32) 

SFR-1 

8.02 

8.03 

36.7 

5.92 

3.04 

(0.  14) 

(0.  17) 

(0.99) 

(0.27) 

(0.40)' 

JRS-1 

7.93 

8.  10 

36.5 

5.77 

2,  17 

(0.  25) 

(0.  30) 

(1.31) 

(0. 43) 

(0.38) 

JRS-2 

8.00 

8.07 

35.  1 

5.63 

3.00 

(0.00) 

(0.25) 

(1.09) 

(0.49) 

(0.37) 

CWS 

8.07 

8.07 

36.7 

5.70 

3.07 

(0.  36) 

(0. 25) 

(0. 90) 

(0.  54) 

(0.45) 

BLR 

3.00 

7.97 

36. 1 

6.00 

3.  13 

(0.00) 

(0.  18) 

(1.07) 

(0. 26) 

(0.35) 

JUR 

8.00 

8.03 

36.5 

5.93 

3.07 

(0.00) 

(0.  32) 

(1.  14) 

(0.25) 

(0.  36) 
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TABLE  8 

VALUES  OF  CHI  SQUARE  FOR  BARTLETT'S  TEST 


Character 

Special  Series 

Santa  Fe  Series 

Main  Series 

D 

30.27*** 

▲ 

67. lO*** 

10.28* 

LL 

21.91*** 

3.69 

SA 

34. 03*** 

28. 56*** 

SB 

17. 04*** 

2*  48 

AG 

2.03 

13.71** 

5.86 

It 

19. 66*** 

37.03*** 

32. 24*** 

21.93*** 

23. 52*** 

58. 14*** 

KL 

mm  mm 

15.71*** 

30.23*** 

35. 94*** 

SL 

20.78*** 

20.80*** 

38. 27*** 

K 

11.88** 

30. 31*** 

31.00*** 

16. 19*** 

35. 78*** 

41. 63*** 

vn 

12. 66** 

23. 66*** 

164. 28*** 

HW 

8.57* 

30.62*** 

46. 25*** 

HL/SL 

54. 32*** 

11.64* 

26.03*** 

HL/E 

16. 23*** 

149.63*** 

37. 02*** 

HL/HW 

16. 79*** 

11. 18* 

37.47*** 

PD/DB 

5.53 

21.49*** 

26. 07*** 

*  P*0.05 
**  P«0.01 
p^  0.001 
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